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Abstract 

We propose a gauged U(1)¥ horizontal symmetry around TeV scale that is a subgroup of a 
SU(3) horizontal gauge symmetry broken at O(10!^) GeV. The breaking generates right-handed 
Majorana neutrino masses through a SU(3)g sextet scalar. A particular Majorana right-handed 
neutrino mass matrix explicitly determines the remnant U(1) 4 at low energy which only couples to 
b—s and u — T in the gauge eigenstate. The dangerous K — K, D — D mixing and B, > utu” are 
kept to be safe because the relevant couplings are suppressed through high powers of small mixing 
angles in the fermion rotation matrix. Our analysis which applies to the general case shows that 
the Tevatron di-muon anomaly can be explained through the B, and B4 mixing while keeping all 


the other experimental constraints within 90 96 C. L. For the B meson decay, the B, + 77 is 


the leading leptonic decay channel which is several orders of magnitude below current experimental 


bound. 


I. INTRODUCTION 


Horizontal gauge symmetry was proposed as an extension of the SM gauge symmetries 
to unify all families of quarks and leptons [1, 2]. Given the three families of quarks and lep- 
tons, SU(3) is the most natural choice for the horizontal gauge symmetry. Interestingly, if 
one assumes all the SM fermions transform under 3 of SU(3)g, the anomaly free condition 
requires three generations of right-handed neutrinos urs [3] while the right-handed neu- 
trinos also play important roles in explaining the origin of neutrino masses. Therefore, the 
SU (3)g horizontal gauge symmetry model provides a natural scheme for the seesaw mecha- 
nism [4] generating small masses for light neutrinos [3]. The Majorana neutrino mass term 
for the right-handed neutrinos explicitly breaks the SU(3)p, thus it is often believed that 
the horizontal gauge symmetry should be broken at a very high-energy scale Mg ~ O(10"*) 
GeV. Then it seems impossible to test the SU(3) gauge interactions in low-energy exper- 
iments. However, it is not always the case as we will show in detail below. Even if some 
subgroup of the SU (3) y remains unbroken, the right-handed neutrinos can still acquire large 
Majorana masses. 

Since ni transform as 3 under SU(3) y, the Majorana neutrino mass term can arise from 
the vacuum expectation value (vev) of an SU (3)p sextet x, 
nig (Xo): (1) 
and Mg = (xs). The light neutrino mass is given by the seesaw mechanism as m, = 
mb(Mg) mp. In order to explain the the neutrino oscillation data, suitable (xg) and mp 
are required. For mp and the other SM fermion masses, there must exist octet Higgs under 
SU (3)g in order to accommodate the correct mass hierarchy in quarks and leptons. In 
addition, to minimize flavor changing effects induced by the octet Higgses, we employ a 
scenario with additional Higgses and singlet fermions [5] in which the mp and quark mass 
matrices or lepton mass matrix are all independent. By taking a suitable gauge choice of 


horizontal symmetry, we always choose the Mg to be a diagonal matrix, 


A00 
Mr=(xe)=|0 B0 (2) 
00C 


The (xs) structure explicitly determines the symmetry breaking. The SU(3) is completely 


2 


broken in a generic vacuum with A Z B # C(# A). However, with a specific vacuum 
of A = B = C for instance, the vacuum (xs) is invariant under a SO(3) symmetry and 
the breaking is SU(3)y — SO(3)y. Being symmetric second rank tensor under SU (3) y, 
the sextet yg transform as yg — UT xgU where U = e'*«T* and T, is the generator of the 
horizontal symmetry. A general scheme to obtain the unbroken symmetry is derived from 
the condition that if (T, (xg) M = 0, (xs) is invariant under transformation defined by T. 
To illustrate the feature of our proposal, we take a vacuum as C = —B. This vacuum 


(xe) = diag( A, B, — B) is invariant under the SU(3) generator Ag as 


000 
T=\5=]001]|. (3) 
010 


Consequently, one can identify the unbroken U(1) gauge symmetry with generator T t and 
it can survive to low energy, for instance O(TeV) which may lead to interesting predictions 
in flavor changing neutral current (FCNC) processes. It was also observed that if there 
exist horizontal gauge interactions, CP violation can be realized with only two generations. 
Explicit examples of CP violation due to U(1)” and SU(2)q was discussed in [1]. If the 
above U(1) is broken at the low energy, the horizontal gauge boson exchanges can induce 
additional CP violations [1] at low energies through quark and lepton mixings. 

In the last decades, huge experimental efforts had been made in improving the measure- 
ments on CP violation in the B meson system. Very recently, the Df) Collaboration at 
Tevatron has reported a large charge asymmetry in like-sign di-muon A^, in both B, and 
Ba decays with 6.1 fb. 
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N N E 


! This U(1) T has an unitary equivalent representation. By taking a 45? rotation R between 2"7 and 374 
axes in horizontal space, the vev becomes 

0 0 

0B]. (4) 

BO 


The U(1) then becomes T" = R'TR = diag(0, 1, —1). We take, throughout this paper, the basis where the 
Majorana mass matrix is diagonal as Mg = diag( A, B, — B) and the unbroken U(1)g generator is given 
by T in Eq. (3). 


where N**(N--) is the event number for bb > u*u* X(u-p- X). Such a large di-muon 
charge asymmetry has a 3.2 ø deviation from the SM prediction At (SM) = (—2.3*0$) x 1074 
[8] and many models have been proposed to account for this anomaly [17-20]. The CDF has 
also measured At, = 8.0 + 9.0 + 6.8 x 107? [21], using 1.6 fb ^! of data, which has a positive 
value and large uncertainties. Combining the above two results in quadrature (include the 


systematic uncertainty), we have 


Ab c —8.5 + 2.8 x 103 . (6) 


sl — 


At the Tevatron both By and B, mesons are produced, hence AP, is related to the charge 


asymmetries ar in B4, and B, decays by ? 
A?, = (0.506 + 0.043)a2, + (0.494 + 0.043)a$, . (7) 


New physics (NP) contributions in B4 mixing are strictly constrained (we will show it more 
explicitly in the parameter fit later), so only large NP contributions to the B, mixing (com- 
paring to the other meson mixings) are allowed. For the NP contribution, if the mixing 
in the rotation matrix between mass eigenstate and gauge eigenstate is not huge, then one 
would naturally expect the U(1)q that maximizes b — s mixing as in Eq. (3). Indeed, for 
a CKM-like rotation matrix, the gauge boson coupling matrix at the tree level in the mass 


basis goes like 


A5 M AÀ 
Gyw] X -X 1 |, (8) 
—A 1 X 


where A is the Wolfenstein parameter [22] around the order of the Cabibbo angle (A œ 0.1). 
Clearly, the meson mixings between the first two generation are highly suppressed. The NP 
also couples to leptons. However, their contributions to the B meson decay branching ratio 


to electron and muon are highly suppressed (although A should be replaced by some small 


2 If the semileptonic b-hadron decays do not involve CP violating phase, then the charge asymmetry is 


directly related to the mixing-induced CP asymmetris in Bg and Bs, meson oscillations. af, = —(1.7 + 
9.1 + 1.5) x 107? [11] is measured through the time dependence of B? + u+ D; X and its CP conjugate 
at DO and ad; = —(4.7 + 4.6) x 107? [10] at the B factory. Nevertheless, we do not use those results here 


due to their large uncertainty. 


mixing of the lepton rotation matrix) ?. Therefore, we focus on the phenomenology in the 
B meson mixing and decay. 

The paper is organized as follows: in section II, we propose the specific model in which we 
consider in the paper. In Section III we show phenomenological implications of our model 
on flavor physics which has subsection III A related to meson mixing and subsection III B 


related to meson decay. Section IV contains our conclusions. 


Il. THE MODEL 


The model starts with a gauged SU(3), model at extremely high energy. By taking all 
the fermions as 3 under SU(3)g. The particle contents under SU(3) jy x SU(2); x U(1)y is 


1 2 
qr : (3,2, 3^ UR: (3, 1, 1), dg : (3,1,—3) 


lr : (3, 2, ESL €R: (3, 1, =2), ng: (3, 1,0) (10) 


which is exactly vectorial and the SU (3) j is therefore anomaly free symmetry. It is crucial 
to have right-handed neutrino triplets, np. for anomaly cancellation [3]. The extension 
to the Pati-Salam unification [23] may be straightforward. 

As we have already discussed the sextet breaking in the introduction, here we focus on 
the Yukawa interactions for the other SM fermions and the Dirac neutrino mass matrix. In 
conventional SU (3); models, in order to break the SU(3) as well as the SU(2);, x U(1), 
one usually introduces one H : (1,2, 1), four ®g : (8,2, 1) to generate all the SM fermion 
mass hierarchies ^. However, the (8,2, 1) Higgs will induce large FCNC [24] if the Higgs 


is light. To avoid the too large FCNC problem, another proposal is to introduce (8, 1,0) 


? Comparing to the other NP coupling matrix in the gauge eigenstate G = diag (1, 1, a) [26], which has a 
coupling matrix 
1 aà — (a—2)33 
G~ aA? 1 (a — 1)? (9) 
(a —2)33 (a— 1)? a 


in the mass eigenstate, the lepton decay branching ratio, especially the one to muon is much more sup- 
pressed in our case. 

4 Another possibility is to consider bulk SU(3) horizontal gauge symmetry broken at one boundary brane 
where all fermion masses are generated at that brane. In this case, one do not need color octet Higgses 
to generate SM fermion hierarchies. See Ref. [27] as an example. 


Higgs|5]. 
i : (8,1,0), H : (1,2,1) “Sade (11) 


In addition, to generate effective Yukawa couplings, a new set of SU (2)z singlet fermions is 


introduced 


4 
2E (3,1, —2), Ex : (8,1, -2), Nx : (3,1,0) 


4 
UR: (3,1, 3); Dg : (3,1, -, Er : (3,1, 22, Ng : (3,1,0) . (12) 


Ur; é (3,1, 


These singlet fermions form invariant Dirac masses and act as messengers to generate the 
necessary Yukawa interactions. We take the up-type quark mass matrix as an example. 
Since the octet Higgs is no longer SU(2); doublet, qrūr® is forbidden and the up-quark 


Yukawa interactions only arise as 
ULO us + MYULUR + qLÜRH + A\ugriirH (13) 


where qruüg.H is universal. After integrating out the heavy fermion fields Uz, Up, the effective 


up-quark Yukawa coupling reduce to 
T(Audiy + (BS) Mg Juh H. (14) 


The same mechanism also applies to the mass generation of down type quarks, charged 
leptons as well as Dirac neutrinos. By assigning the (0) independently, the mixings 
and masses in different fermion sectors are completely independent for each other and one 
can easily accommodate hierarchies and mixings in SM fermions and the Dirac neutrinos. 
This also enables us to choose the Dirac neutrino mass matrix other than nearly-diagonal 
structure. 

After electroweak symmetry breaking, the effective Yukawa coupling of ow also arises 
as 


(H) My tipo u; . (15) 


Then, the oO exchanges induce FCNC’s in general. We have checked that they satisfy the 
strongest constraint from K-K mixing, marginally ?. However, actual effects depends on the 


mass spectrum of the oO and hence we do not discuss them in this paper. 


5 Since M ~ (®g) = Mz://gg ~ 50 TeV, the effective coupling here is (H)/M of O(107?). With additional 
propagator suppression due to 1/M?, the amplitude is O(10- 19) GeV ?. 
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Another consequence is that both up and down quark mass matrices become Hermitian 
mi = m, mj, = ma (16) 


Thus, the CP violation in strong interactions due to quark mass matrices, 
arg(det(m,) det(ma)} is absent at least at the tree level[5]. In addition, the Hermit mass 
matrices also require the rotations Ur, Ug in the mass diagonalization Ul m,U. r to be equal 
Uz = Ug. In this case, the horizontal gauge boson couples to vector currents of quarks and 
leptons. As a consequence, pseudo-scalar bosons like B, or Ba do not decay to a pair of 
leptons. However, this is only the result of our specific choice of mass generation model for 
quarks and leptons. In the following analyses, we assume more generic rotation matrices 


and Ug and U; are taken independent for each other to estimate the prediction. 


III. PHENOMENOLOGICAL IMPLICATIONS IN FLAVOR PHYSICS 


The horizontal gauge interaction is real but family dependent. After the mass diago- 
nalization, the other flavor violation entries as well as new CP violation can arise. The 


Lagrangian of gauge interactions is 


JE Lu = gu LT" Za + LOR 
gud, (V'TVI) Hah Z +L eR, (17) 
ij 


where V7? stands for the rotation for left-handed q-type quarks and T is the generator of 
U (1)g interaction given in Eq. (3). 

Flavor changing interactions in the SM can only be measured via electroweak charged 
current interactions. Therefore, for the SM fermion rotation matrixes, only the left-handed 
ones get constrained from the CKM matrix V“(V?)' = Vex and one cannot determine 
even V? and V£ respectively. The other rotations are completely unknown. For simplicity 
of the discussion here, we will assume that all magnitudes of the left-handed mixings are 
CKM-like but the complex phases are O(1) and unconstrained right-handed mixings have 


the similar structure. Therefore, we have the mixing matrix in the mass eigenstates as 


Am. De. LEX 
(Gs = VEIT (VE) a G e| o 7 1 (18) 
eA d X 


This U(1) gauge interaction maximizes the mixing in between second and third generations. 
Mixing magnitude in B,, By and K” is at the order of (1: A? : A). The D? — D? mixing 
is also at the order of A9 suppression comparing with B, mixing. This U(1) is consistent 
with the phenomenological constraints among different meson mixings. If one assume the 
lepton doublet and right-handed singlet rotations are the similar to the quark sector °, one 
can also compute the leptonic decay of mesons. For instance, B, > utp decay partial 


width has a A* suppression. 


A. Meson Mixing 


At the energy scale mp, the effective Hamiltonian responsible for neutral meson mixing 


(and in particular B, — B, mixing) through the tree-level exchange of Z' is 
H= Cj, (m) OF, t Cal) ORR + Cpl) O7n T CYR) p (19) 
where the AF = 2 operators are given by 


OF, = dry" Prajd/y, Pra; Og = di" PrljliYu PRI, 


Ofr = GY" Pig GpPry, Ofr = &PLyGPry. , (20) 


and the Wilson coefficients at Mz scale are (C9. (M. z') = 0) 


2 2 2 
ij g lij ij g lij ij g tij 
Ce, (Mz) = (GP) COhr(Mz) = E (GYP C#,(Mz) = (GP) (21) 
M2, M2, 


where gy is the horizontal gauge coupling at Mz, scale and Mz, is the horizontal gauge 
boson mass. 

In order to calculate the B physics observables, one has to take into account the running 
effects of the four operators above. The relation between these four operators at the Mz 


and m, scale is presented in Appendix A. After one obtains the Wilson coefficients at my 


$ Within the minimal SU (5) or SO(10) grand unification theory (GUT), both left-handed and the conjugate 
of right-handed states are embedded into the same GUT multiplet. But in the SU (3) x, left-handed states 
and right-handed states both transform as 3. The GUT multiplet will contain both 3 and 3 and the 
horizontal gauge symmetry model is not consistent with the minimal GUT. Therefore, we will not assume 


any correlation among rotations of quarks and leptons. 


scale, by using the relevant hadronic matrix elements [39] 


= 1 
b b 
(B4|O 7, rel Pa) 3B, f, BiL/RR 
= 1 
(BylOpa|Ba) F -zme fB, BiR 
"IE 1 " 
(BOW |B) = cma, 3, BU (22) 


2 Dens bd x Bt ~ pbr — ple — 
Here we use m5 /(m, + mq)” © 1 and assume By, ~ Bir ~ Brr = Bp, = Be, Then we 


obtain Mh 


- 1 C" (u) 369 (y 
Mf, = (Bal HIB,) = — 5 f2,ma, Bo, | CI Q0) + Clu) - Ces) 4 SCR UU (us 


From the discussion above, the flavor off diagonal coupling between the horizontal gauge 
boson Z' and the first two generation quarks are highly suppressed (For the CKM like 
rotation matrix, it is at least A? suppressed), so we will neglect the new physics contributions 
to the K — K " and D — D mixing. The Z'-b-s and Z'-b-d couplings, on the other hand, 
is either unsuppressed or A suppressed, hence we expect large new physics contributions to 
modify the magnitudes and the phases of M1/^, where M17? are off-diagonal mixing matrix 


elements in Eq.23. We can parametrize such effects by ? 
s s io^ 
M = (MG) Aas As = (Aale. (24) 


The experimental measured observables are summarized as follows, Ama, and AT a/s 


measures the mass and decay width difference between the heavy and light mass eigenstates 


9 „d/s 


of the Bays mesons ”. ag; is the charge asymmetry in semileptonic Bajs decays. Ba or B, 


measure the time-dependent CP violating phases in the hadronic B decay channel By — 


T The constraint from the CP violation in the K — K mixing [28] is, in fact, marginal for A 0.1, gg œ 0.02 
and Mz; ~ 1 TeV [I1]. 

8 The CP violation phase is defined as ¢ = Arg(— Mi2/T42) so we choose T, to be real here. 

? The experimental uncertainty in the measurements AT'4 = sign(ReAcp)AT'4/T4 = 0.009 + 0.037. is too 
big for us to consider it. 
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FIG. 1: The Abs(G?.) versus Arg(G7.) region plot for Mz = 1 TeV and gg = 0.02. In the left 
panel, we choose Abs(G,;)/Abs(G;.) = 0.1 and Arg(G;,) = Arg(G;,) while in the right panel, 
we choose Abs(G;,)/ Abs(G;.) = 0.05 and Arg(G;,;) = Arg(G;.). The magenta, green, and yellow 
region with dashed, dot-dashed and dotted boundary stands for the allowed parameter space at 
9076 C. L. for Amsja, ATs and Syx/s. 
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FIG. 2: The similar plots as FIG. 1 but combine both the B4 and B, experimental constraints. 
The small overlapped region labeled with “All” is the parameter space that fits for all experimental 


constraints at 9096 C.L. 
J/v Ks or B, > J/v à. They are shifted by the CPV phases in Ba or D, mixing. 


Ama/s = Am‘) |Aays| » 


AT, = AT™ cai" ado). 
d/s AT sin(o5)5 + Pijs) 


a —4 
SL = Ams Aga 
Sux = sin (267+ 4), 
2prP = 2b- 92. (25) 


The theoretical inputs are listed in the Table I. All the decay constants and bag param- 
eters are used from Ref. [37]. Notice that we use the calculations for Amp, = 2|M?,|°™, 
Amp, = 2|Mj,|5“ and $? in Ref. [36] which uses the more recent decay constants and bag 
parameters with much smaller uncertainties. All of the rest SM inputs are either from [8] 
or [6] 

All experimental measurements which are used to compare with our model outputs are 
listed in Table II. For the like-sign dimuon charge asymmetry, we use A^; ~ — (8.5 + 2.8) x 
107? which combine the DØ measurements with the CDF measurements. For GF? and T, 


measured by both CDF and DØ [12-15], we use the combined results with each measurements 
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TABLE I: The theoretical input parameters [8]. 


Ba Bs MBa MB, 
0.38 + 0.01 [9] 0.018 + 0.001 5.28 GeV 5.37 GeV 
fp, V Ba fs, V Bs fBa fB. 
(216 + 15) MeV (275 +13) MeV 192.8+9.9 MeV 238.8 + 9.5 MeV 
oy oe (Amp,)™ (Amg, )™ 
—0.091*0 038 (4.2 + 1.4) x 107? 0.53 + 0.12 ps! 19.30 + 2.2 ps! 
(AT4)P*P (Are? 
(2:07 7025 x 1073 ps7! 0.098 + 0.024 ps! 


in 2.8 fb! [14] where do not add the most recent CDF results [15] since it is difficult for us 
to extrapolated their contributions between 2.8 — 5.2 fb^!. Nevertheless, we find that our 
conclusion in the paper does not change if we use the combined results in Ref. [20] 1°. 

The experimental constraints on the parameter space of our model are presented in Fig 1 
and 2. The parameter Abs(G;.) and Arg(G;.) are quite similar as the parameter h? and c, in 
Ref. [25] except for a overall factor with very small phase related to the (M%,)°™, therefore 
our results here can be used as the model independent analysis after some re-parametrization. 
For illustration, we choose parameter for By couplings Abs(G5;) = Abs(G4,,)/10, Arg(Gi,4) = 
Arg(G^.) in which there are sizable contribution to A^, from B4 and Abs(G‘,,) = Abs(G;.)/20, 
Arg(G,;) = Arg(G;.) in which the af, contribution to A^, is negligible. In contrast to the 
paper [18], it is clearly that from the upper right plot in FIG 1, there is no region allowed 
by all experimental constraints within the 1 o !'. The best fitted region for the phase 
Arg(G,.) C (1/2,37/4) are quite consistent with the one found in Ref. [25]. However, since 
the Ref. [25] essentially marginalize over T}? in the range 0 — 0.25 ps^' and use the best fit 


points in which AT, is about 2.5 times larger than the prediction, the goodness of fit in our 


10 For the BE*P at 9096 C. L., we find BE*P c (0.13,0.74) from Fig. 7 in Ref. [14] since 1 - CL is almost 
$ g 


linear to Log(L) where L is the likelihood ratio. 

11 We notice that in [18], they use SM prediction directly from [8] which use the old decay constants and bag 
parameters. Nevertheless, we find that our conclusion still holds because the parameter space mentioned 
in [18] are not allowed by (Amp,)P*P and (AT,)E*P. In the other paper which contains the model 
independent fit [19], they simply neglect all the large uncertainties from the SM prediction which lead to 
the wrong upper bound of A}. 
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result here is reduced significantly comparing to the one in Ref. [25]. 


TABLE II: The experimental data. 


STK (Bs) EP (Amp,)*P (Amp, P"? (Ar)? 
0.655 + 0.024 [10] 0.441017 0.507 + 0.005 ps! 17.77 +0.12 ps™! 0.154+9:054 pg 7 
Ab BE*P(B, — ete—) BE*P(B, > utu-) BE*P(B, — et pt) BP*P( B0 — prt) 
—(8522.8 x10? ps! «54x10 < 4.7 x 1078 < 6.1 x 10-9 < 3.8 x 1075 


B. Meson Decays 


The new horizontal gauge boson can also mediate meson hadronic decay or leptonic 


decays at tree level. In this session, we choose to discuss the two leading processes, b — 


scc and b + spur respectively to illustrate how meson decays constrain the horizontal 
gauge interaction. The other transitions are always with additional factors of power of A 
suppressions. 


The effective AF = 1 Hamiltonian responsible for neutral meson decay is 
HAT = C3Q3 + CsQs + C3Q3 + CsQs (26) 
and the leptonic decay like Bg, > ¢T¢-, 
Hig = CoQo + CioQio + CoQo + CioQio (27) 


where 7, j = e, pu, T as lepton flavor indices. 

In the case of b — scc transition, SM contribution at tree level is induced via weak 
charged current with a CKM factor V4,V*. ~ A?. Reading from the effective couplings, 
the horizontal gauge boson mediated b — scc has a factor of GysGee ~ A?. The SM and 
horizontal gauge interaction contributions are at the same order of A and the one can simply 
compare their couplings and gauge boson masses to estimate the ratio. As we discuss in 
previous section, the U(1)g is broken at O(TeV) which results in a suppression due to 
(gu /g)*(miy/ M) ~ 1075. Consequently, the contribution to b > scc from new horizontal 


gauge boson is completely negligible. In the case of neutral meson mixing, SM leading 


contribution is from box-diagram while the horizontal gauge boson contribution is at tree 
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level. Therefore, even if the new horizontal gauge boson is of order TeV, it is still possible 
to change the AM significantly. For decay process, if there exists SM tree level, the above 
argument then always applies. We won’t discuss any constraint from such decays. 

Since the horizontal gauge boson also couples to the leptons, there is again tree level 
contribution to the meson leptonic decay. Within the framework of SM, B, pure leptonic 
decays are realized via the electroweak penguin diagrams with Z/4* — (*(^ and therefore, 
there is no lepton flavor violation at all. The constraints on leptonic decay are mostly on 


leptons directly decaying from B meson. If there exists a 7 in the final state, 7 decay will 


complicate the search due to the D* decays. Therefore, the leading constraints are 


Br(B,— utu) < 4.7 x 107° 


Br(B, 2 e*e^) < 5.4 x 107° 


Br(B, => e*u?) < 6.1 x 107° (28) 


In our model, the horizontal gauge boson has maximal coupling to b, s and u,T and the 


leading leptonic decay constraint is from B, — u-*T*. But as we mentioned, B, — u^r* 


does not exist in SM physics at leading order and it is only from the new physics contribution. 


Using(0 | sy"“75b | B9) = «fpp'5, one can compute the decay BR as 
B 


3 2\ 3 2 
Br(B, > u^r*) = Tal (B° — yt) = fs, Tor ¢ — i) (=) Chia (29) 


MB 
where the Wilson coefficients Co19 = g3(m»)/m%,. To estimate the decay BR, we take 


Mz ~ 10? GeV, gy œ 0.02 and substitute m, = 4.7 GeV, m, = 1.777 GeV, fp, = 230 MeV, 
mp = 5.3 GeV, Tg = 1.6 ps, we obtain 


Br( By => u*rF) ~ 1.2 x 107°. 30 
H 


12 The only relevant search is from CLEO as of B4 > u*7* [30], 


Br(B4 — u*7*) < 3.8 x 107°. (31) 


Due to horizontal gauge boson given, the By decay partial width has an additional A? 


suppression so the bound is prediction well below the experimental bound. One can also 


12 The result is based on assumption from left-handed chiral interaction only which is sufficient to estimate 


the maximal value of leptonic decay BR. 
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estimate the B, — utu using the above result. The B, — pty has a factor of M 
suppression then the prediction is about two orders lower than the current experimental 
bound. 


The other possible rare decay which can be induced by the horizontal gauge boson is the 


FCNC decay in top quark, for instance, t > c/u + u^7r*. However, given the large Mz, the 
three body decay is highly suppressed. 


GLm) (my \* f gg M. 
T(t > c/u + utr*) = F £3 ( 3 expe 103 (32) 


Even at the top factory like Large Hadron Collider, it is impossible to observe such rare 


decay event. 


IV. CONCLUSION 


The dimuon asymmetry reported by D@ Collaboration is much larger than the SM predic- 
tion which suggests new sources for CP violation. In this paper, we propose the possibility 
to explain such an anomaly through a tree level exchange of a gauged U(1)g horizontal 
symmetry in B meson mixing. The U(1)y horizontal symmetry is a remnant symmetry of 
SU(3)g broken at Mg ~ O(10™) GeV through a sextet scalar which gives the neutrino 
mass. Such a U(1)g gauge boson only couples to b — s and u — T in the gauge eigenstate 
which suppresses all other dangerous meson mixings and B meson decays after the flavor 
rotation. For a general flavor rotation matrix we consider there is a parameter region around 


the phase Arg(G;.) C (7/2,3-/4) which fits the data at 9096 C. L. For the B decay, the 
bs 


dominate enhanced channel is the B, — u^7*. Nevertheless, such a enhanced decay channel 


is still one order of magnitude smaller than the current experimental bound. 
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Appendix A: RG Running of the AB = 2 Operators 


In the previous section, the Wilson coefficients are given at the Ms scale while to calculate 
the physics processes involving low energy mesons, one will need to calculate relevant Wilson 
coefficients at the low energy scale. The running contains two steps, the first step is from 
Mz scale to m, where six flavors contribute to the running of a,, the second step is from m; 
to mp, where only five flavors contribute. We summarize the running effects of the relevant 
AF — 2 operators below [39]. 

The operators belonging to the LL/RR, LR sectors read 


Qu, = (s*y, Prb) (s? PP), 
Qrr = (so, PLU?) (s^ Prb), 
Qua = (8° I Pub?) ("4^ Prb"). 
Qua = (s*PLU*) (s? Prb). (33) 


Cur(uo) = [n(Ho)|Crr (Hr), 


Cir(Hs) ma(ts)] [maus] Crn(u) 


4 du (34) 
Crn(u) nai(Hs)| [ros (uo)] Crn(u) 


( 

Crr(Hs) = [n(uo)|Cnn(u), 
| 
| 


The variables 7; is defined as ratio between strong coupling constant at different scales. 
Given the U(1) gauge boson is around 10 TeV, we have two different running of a, taking 
into account the threshold correction due to top quark. For the running between the B 


physics scale and top quark, we have 


(35) 


For the explicit form of 7, we list both the expression at the LO (with subscript (0)) and 
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NLO (with subscript (1)) 


Similarly, for the 


boson scale, one can 


nO iy) as 

n (u) = 1.6273(1 — ns)ns , 

(HO ce aa 

ni? (ue) = 0, 

na) = s(n ago 

vo (ub) coeno s 

n (uy) = 0.9250; 24/3 + 3/2 (2.0994 + 1.17445) , 

mp (He) = 1.8875(q5 7 — g 24/22) , 

nSP (u) = (—11.7329 + 0.7829;)n2/?? — ng“ (—5.3048 + 16.25485) , 

nS (ts) = (7.9572 — 8.8822ns n; ^? + 0.925020 . (36) 


running between the top quark mass scale and the horizontal gauge 


replace all the py, pi by pe, HM, in Eq. (34), where all the 7s are 


a tap) eme s 
nO (uy) = 1.3707(1 = ni , 
nuo eqs 
"i2 (up) — 0, 
P(e) = sre = a), 
ee (ie) Se 
nU (uy) = 0.921974?! + ng (2.2194 + 1.29759) , 
Jha (i5). — 98980 eqs 7) 5 
nU (uy) = (—10.1463 + 0.865095) ne + ng“ 6.4603 + 15.7415) , 
nS) (uy) = (9.6904 — 10.612255); ^!/?! + 0.921952?! . (37) 
and a (u) 
n = 6 (38) 
os (uu) 
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